Time-independent approximations for periodically driven systems with friction 
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The classical dynamics of a particle that is driven by a rapidly oscillating potential (with frequency 
Lo) is studied. The motion is separated into a slow part and a fast part that oscillates around the 
slow part. The motion of the slow part is found to be described by a time-independent equation 
that is derived as an expansion in orders of uj~^ (in this paper terms to the order lu~^ are calculated 
explicitly). This time-independent equation is used to calculate the attracting fixed points and their 
basins of attraction. The results are found to be in excellent agreement with numerical solutions of 
the original time-dependent problem. 
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Time-dependent systems typically exhibit behavior 
which is more complicated than the one of the corre- 
sponding time- independent ones. Moreover, physicists 
did not develop yet an intuition on the dynamics of time- 
dependent systems to the level of the one that exists for 
time-independent ones. As a result, it is difficult to pre- 
dict the qualitative properties of driven systems even in 
cases in which it is easy to understand the dynamics of 
time- independent systems that are similar. In this work, 
the dynamics of some time-dependent driven systems will 
be related to the dynamics of time independent ones. 
In particular, this will be done for some systems where 
there is a clear separation of time scales. This will enable 
qualitative and quantitative analysis of the dynamics of 
some class of time-dependent systems using the experi- 
ence with time-independent ones. 

There is a vast literature on methods which approx- 
imate a dynamical system by a "simpler" system with 
smaller number of degrees of freedom. These include, 
among others, averaging methods 1, 2], multiple time 
scale analysis j^i, and centre manifold theory |^]. In fact, 
the method employed in this article can be viewed as a 
particular example for the more general method of mul- 
tiple time scales analysis 2] . It should be noted that the 
method used in this work is tailored for equations of the 
form (Q, and therefore its application is much simpler 
then the use one of the more general methods. However, 
it is representative of many physical problems. 

The more known case of separation of time scales is 
the adiabatic one. In this case, the system evolves on a 
time scale which is much shorter then the time scale of 
the driving. There are many works which treat different 
adiabatic approximations 0, 0- On a qualitative level, 
one can understand the resulting dynamics by treating 
the driving as if it was fixed in time, letting the system 
evolve, and then change the values of the parameters of 
the driving according to their time dependence. 

Less known, but not less interesting, is the opposite 
limit. In this case, the typical time scales of the dy- 
namics of the system in absence of the driving are much 
longer than the period of the driving. A remarkable ef- 
fect of such driving is "dynamical stabilization" , in which 
a particle, that in absence of the driving can escape from 



some region, may be trapped by the rapidly oscillating 
field. Exarnples for this phenomenon include the Kapitza 
pendulum and the Paul trap 

A simple treatment of rapidly driven systems is given 
in textbooks j^] , following Kapitza's work on the inverted 
pendulum 0. In this approximate calculation the mo- 
tion is separated into a sum of a slow part and rapid 
oscillations around it. The rapid oscillations are com- 
puted explicitly and their effect on the slow motion is 
found. The treatment of Landau and Lifshitz turns out 
to be the leading order of an expansion in (where u) 
is the driving frequency)- It was extended to the order 
Lu^^ in recent works In these works it was demon- 

strated that for rapidly driven Hamiltonian systems it is 
possible to obtain a time independent Hamiltonian that 
controls the slow motion. It is obtained by a canonical 
transformation as an expansion in powers of oj~^. 

In this article, the method developed in 0, 13 is ex- 
tended to rapidly driven classical systems in the presence 
of friction where the Hamiltonian formalism is inappli- 
cable, since energy is dissipated. The (modest) goal of 
this article is to demonstrate how one can understand 
and predict the qualitative properties of the dynamics 
using very simple approximations. In particular, it will 
be demonstrated how this method can be used in order 
to predict qualitatively the general form of the basins of 
attraction for such systems and how to compute their 
boundaries from the equation of motion of the slow part. 

Newton's equation of motion, for a rapidly driven sys- 
tem, is given by 



mx + ax — —Vq{x) — V{{x, ujt), 



(1) 



where m is the mass of the particle and a is the friction 
constant. The potential is y(x,wi) = Vo{x) + Vi{x,ujt) 
that is chosen so that the time average of Vi over a period 
vanishes. Derivatives with respect to the coordinate and 
time oi f{x, t) are denoted by /' and / respectively, while 
/ denotes the average over a period. 

The slow and fast motion are separated with the help 
of the ansatz 



x{t) ^ X{t)+i{x,X 



(2) 
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where ^ is periodic in the fast time variable r = ujt, with 
a vanishing average over a period, and X{t) is the slow 
part of the motion. The fast time r is treated as an 
independent variable. One ensures that X{t) is indeed 
slow by choosing ^ in such a way so that the equation of 
motion for X does not depend explicitly on time. This 
can be done, at least approximately, by expanding ^ in a 
power series in uj~^, using 

n=l 

Then, the functions ^„ are determined, order by order (in 
0J~^), from the condition that the remaining equation of 
motion for X is time-independent. This procedure leads 
to new time-independent terms in the equation of motion 
for X{t). (These terms cannot be canceled by a choice 
of which stay bounded for large times.) An explicit 
derivation of the first few terms for a similar problem 
can be found in [^. For an equation of motion of the 
form JQ) , the perturbation theory results in the following 
explicit solution for ^ (given in terms of X{t)) 

/ . \ 1 r^^^^ 2X 

\ / muj^ J muj'^ J 

+ / [V;]+0{u;-'). (4) 

Substitution of Q in Q leads to the slow equation of 
motion, for X{t), 

mX = -aX - V^iX) -^f [¥{'] f [¥{] 

mw^ J J 

+ ^^/[^"J f'\yi\+o{u.-% (5) 

The symbol [...] denotes an integral over r, defined 
only for periodic functions of r with vanishing average, 
which is performed in such a way so that the resulting 
function of r is also periodic with vanishing average. This 
integral is easy to compute using the Fourier expansion of 
the integrand (see 131 for details). Multiple apphcation 
of the integral (j times) is denoted by J^"'"'^ [...]. The 
overline denotes an average over a period of r. There- 
fore, equation (|SJ| does not depend explicitly on time (and 
its solution will not exhibit oscillations with the external 
frequency uj). 

The leading-order (w~^) correction to the motion in 
absence of the driving can be seen as resulting from the 
effective potential 

K//(^) = V,{X) + ^ (^j^ [Vi{X, r)]) (6) 

that acts on the slow motion. However, at higher or- 
ders, terms appear which do not seem to result from a 
potential. Note that in spite of the rapid oscillations, the 



friction, and hence the energy dissipation, is associated 
only with the slow motion (at the order lo^"^). 

Equations Q and (0) are the results of a high fre- 
quency perturbation theory. They result in a mapping of 
a time-dependent problem into a time-independent one. 
In any given order in uj^^, this theory reproduces the 
result of the mathematical theory of separation of time 
scales, but it is much simpler. One of the goals of this 
article is to demonstrate how these equations can be used 
to understand the dynamics of such driven systems. It 
will be demonstrated by a simple example. 

Consider a particle which, apart from the friction, is 
under the influence of an oscillating field given by 

Vi [x, ut) = Ae"^^'' sin(fca;) cos(t^i -f 0). (7) 

This simple system is of interest since the time average 
of the potential vanishes, Vq{x) = 0. Therefore, the in- 
fluence of the oscillatory field is dominant even at high 
frequencies, in contrast to systems with Vq{x) ^ 0. Ac- 
cording to the high frequency perturbation theory, the 
slow part of the motion of the particle can be (approxi- 
mately) viewed as motion, with friction, in the effective 
potential 

Veff{X) = [/ccos(fcX) - 2/3Xsin(fcX)]^ . 

(8) 

This effective potential is depicted in Fig. ^ It exhibits 

0.012 I ' , ' , ^ , ^ , ^ , ^ , ^ , ^ 1 
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FIG. 1: The effective potential for the slow motion. The 
values of the parameters in Q are A = m = l,/3 = 4, fc = 2 
and cj = 10. 

several minima. These are of interest since a particle 
moving in a time-independent potential, in the presence 
of friction, will be found at one of these minima after a 
sufficiently long time, in spite of the fact that the linear 
stability of these points is time-dependent in the original 
system. While the potential JHl, depicted in Fig. ^ has 
several minima, two of those, at x ~ ±0.3266, are more 
pronounced. Based on these observations, one can pre- 
dict that after sufficiently long time the particle will be 
found in the vicinity of one of the minima also for the 
time-dependent system. 
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FIG. 2: Basins of attraction of x(t oo) ~ ±0.3266 for the 
time dependent system Q. Results for several frequencies 
are presented, while A = m = l, l3 = 4, k = 2 and 4> = 0. 



In this article, we are interested mainly in the basins 
of attraction of the two minima located at x ~ ±0.3266, 
that is, the initial values xq = x{0),vo = i(0) which 
evolve to these minima for long times. These basins of 
attraction, for different values of the frequency, are de- 
picted in Fig. 121 The basins of attraction exhibit some 
qualitative properties which are of interest. At lower fre- 
quencies {uj = 6, 10) these basins are separated while 
for higher frequencies (w = 20, 40) they have a common 
boundary. This qualitative behavior can be understood 
with the help of the effective potential (Fig. ^| for the 
approximate time-independent system. 

Consider the trajectory (of the effective time- 
independent system) which starts with the initial val- 
ues Xq = 0, Xq — £ > (that is, arbitrarily close to 
Xo = Xq = 0). If this trajectory, after a long time, 
manages to pass over the first maxima of Vef / with pos- 
itive X, then any other trajectory which starts say with 
Xq < and a (large enough) positive velocity will also do 
so. In this case, the basins of attraction of x ~ ±0.3266 
are separated. In contrast, if this orbit, starting at 
Xq = and Xq — e > 0, is trapped at x ~ 0.3266 after a 
long time then the basins do have a common boundary 
which includes the point Xq = Xq = 0. Consider the 
effective potential depicted in Fig. ^ A particle located 
near X = will feel an effective force due to it. This ef- 
fective force will accelerate the particle while the friction 
will decelerate it. If the friction is dominant the particle 
will be trapped in the first minimum and it will be found 
near x — 0.3266... after a long time. In contrast, if the 
effective force is strong enough one can expect that the 
particle will be able to pass the first maximum of Vgf / (X) 
at positive X and therefore will not be near x ~ 0.3266 
after a long time. In the example used in the present nu- 
merical investigation the friction does not scale with the 



frequency while the effective potential (and the resulting 
force) scales as Therefore, one can expect to find a 
transition from separated basins of attraction to basins 
with a common boundary as the frequency is increased. 
This is consistent with the numerical results presented in 
Fig. El 

It was demonstrated that some qualitative properties 
of the time-dependent system are described by properties 
of an appropriate approximate time-independent system. 
It is of interest to see whether one can obtain also quanti- 
tative results using the high-frequency perturbation the- 
ory. The boundaries of one of the basins of attraction 
of a; ~ 0.3266 of the time-dependent system are com- 
pared to the ones resulting from the approximated time- 
independent system in Fig. 13 The motion of the time 
dependent system depends also on the phase of the os- 
cillating force at i = 0. It is important to note that this 
comparison is fairly naive since the initial values xq and 
vq are only approximately equal to Xq and Xq of the 
time independent system. Fig. |31 shows that the bound- 
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FIG. 3: The boundary of the basin of attraction of a: ~ 0.3266 
of the time-dependent system, for various phases, (thin lines) 
are compared to the one of the time-independent effective 
system (heavy line), for different frequencies. The values of 
the phase are = 0, n/2, tt and 3tv/2. 

ary of the basin of attraction of the time-dependent sys- 
tem fluctuates around the corresponding boundary of the 
effective time-independent system when the phase of the 
field is varied. This is not surprising since it is clear that 
the basins of attraction must depend on this phase. 

As was mentioned earlier, the comparison in Fig. |3| 
is naive, in the sense that the coordinates x{0),v{0) are 
not the same as the coordinates X{0),X{Q) of the time- 
independent system. However, the results presented in 
Fig. O still demonstrate that at high frequencies, the 
basins of attraction of the time-dependent system can 
be approximated by those of time-independent effective 
ones. The size of the fluctuations seem to decrease when 
the frequency is increased. 
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To obtain a better quantitative correspondence be- 
tween the results obtained using the high-frequency per- 
turbation theory and the numerical results for the time- 
dependent system, one has to account for the difference 
between the slow coordinates X,X and x,v. The con- 
nection between x{t) and X{t) (and also X{t)) is given 
by Eq. Q • One can use the high frequency perturbation 
theory to obtain an expansion for ^ and then substitute 
t — 0. This leads to an equation for the initial value 
xo = x{0) in terms of Xq = X{0),Xo = X{0). Simi- 
larly, by differentiating with respect to time at < = one 
obtains vq = v{0). The calculation results in 



merical results, presented in Figs. |51 and |21 demonstrate 
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e-P^'o cos(0)/(Xo) + 0{uj-^) 



Xo e-f^'o sm{<j>)f{Xo) 



aA 



e-^-^oXoCOs(0).g(Xo) 

-'3^0 cos(0)/(Xo) + O{uo--') 



(9) 



where /{Xq) = kcos{kXo) — 2/3Xo sin(A:Xo), while 
g{Xo) = [4/32 -P- 2(3] sm{kXo)-4:(3Xokcos{kXo). 
From Eq. Q it is clear that the fluctuations in vq are 
of order uj~^ when the phase is varied, while the fluctua- 
tions in Xo scale only as a;"^. This is in agreement with 
the fluctuations presented in Fig. 13 

To test the perturbation theory more quantitatively, 
the boundary of the basin of attraction of the effective 
time independent system was mapped back to the orig- 
inal coordinates using Eq. @ and compared to the nu- 
merical results obtained for the time-dependent system. 
To avoid complicated graphs, we only present the com- 
parison for the phase = tt/2. The results are presented 
in Fig. ^ It is clear that the agreement is excellent. 
Only a small difference is seen for lo = 6 while for w = 20 
the difference between the boundaries cannot be observed 
on the plot. This demonstrates that the high-frequency 
perturbation theory can be used to obtain quantitative 
results and not only qualitative ones. 

In this article, we have used a high frequency pertur- 
bation theory to describe the motion of a rapidly driven 
classical particle in the presence of friction. In this per- 
turbation theory the motion is separated into a slow part 
and rapid oscillations around it (see Eqs. 0-©). An 
equation of motion for the slow coordinate, accurate to 
order u;^'^, was obtained. The oscillations of the particle 
around the slow solution were calculated as well. The 
slow motion is found to be approximately described, to 
the order w"^, by the motion of a particle in an effective 
potential with friction. This suggests that after a long 
time the particle will be found at the minimum of this 
potential, a fact which is confirmed numerically also for 
the time-dependent system given by Eq. {Tj) . Due to the 
dissipation, these minima are surrounded by basins of 
attraction, which include the initial phase space points 
that flow to those minima after a long time. The nu- 




FIG. 4: Comparison between the boundary of the basin of at- 
traction of s ~ 0.3266 of the time-independent system (heavy 
line) and the boundary of the time-dependent one (thin line), 
for different frequencies, where (j) — 7r/2 while the rest of the 
parameters are those of Fig. |5| 



that some of the qualitative features of the basins of at- 
traction can be understood by considering the simpler 
time-independent effective system with the potential of 
Fig. n It was also shown, in Fig. 0] that, by a more 
careful analysis, one can obtain excellent quantitative 
agreement between the motion generated by the time- 
dependent equation of motion and the motion generated 
by the corresponding time-independent one. 

The results presented in this article suggest that the 
high frequency perturbation theory can be used to ob- 
tain both qualitative and quantitative understanding of 
the dynamics of a classical rapidly driven particle in the 
presence of friction. The friction effectively dissipates the 
energy (only) of the slow motion (up to the order 
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